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ABSTRACT 

X-ray reflection spectra are an important component in the X-ray spectra of many 
active galactic nuclei and Galactic black hole candidates. It is likely that reflection 
takes place from highly ionized surfaces of the accretion disc in some cases. This can 
lead to strong Comptonization of the emergent iron, and other, absorption and emis- 
sion features. We present such reflection spectra here, computed in a self-consistent 
manner with the method described by Ross and Fabian. In particular we emphasise 
the range where the ionization parameter (the flux to density ratio) ^ is around and 
above lO'^. Such spectra may be relevant to the observed spectral features found in 
black hole candidates such as Cygnus X-1 in the low/hard state. 
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1 INTRODUCTION 

Compton reflection is an important component of the X-ray 
spectrum of many compact objects (Guilbert & Rees 1988; 
Lightman & White 1989), particularly of Active Galactic 
Nuclei (AGN) and Black Hole Candidate (BHC) sources. 
It is produced when the primary X-ray emission from the 
source strikes Thomson-thick matter which Compton scat- 
ters X-rays into our line of sight. In many cases the accretion 
flow, probably a disc, is the major scattering medium. If that 
matter were solely composed of hydrogen then the reflection 
continuum would have the same spectral shape as the pri- 
mary emission in the X-ray band, decreasing above about 
30 keV due to Compton recoil. Other elements such as oxy- 
gen and iron, if they are not highly ionized, can absorb the 
lower energy X-rays and so flatten the reflection continuum 
and, in particular iron, add fluorescent line emission (George 
& Fabian 1991; Matt, Perola & Piro 1991). Measurements 
of the strength and shape of the reflection spectrum can 
yield the geometry, velocity, gravitational potential depth 
and abundances of the scattering medium. 

As the X-ray irradiation of the medium becomes more 
intense, so the matter becomes ionized. This reduces the 
effects of absorption, progressively from lower to higher en- 
ergies as the intensity increases (Lightman & White 1988; 
Done et al 1992; Ross & Fabian 1993; Zycki et al 1994). 
An important signature of ionized reflection is a strong iron 
edge. Iron is a strong absorber when the gas is neutral, but 
absorption both above and below the edge are so strong that 
the change in observed flux across the edge is small, par- 
ticularly if the primary continuum is also observed. When 
much of the oxygen in the gas is completely ionized the 



edge appears much stronger because the absorption below 
it is reduced. Of course the energy of the iron emission line 
shifts up in energy from 6.4 through 6.7 to 6.97 keV as 
iron is ionized, but doppler shifts may confuse precise es- 
timates of the observed energy. The line can be weak when 
iron is ionized to Fe xvil-Fe xxill since line photons can be 
resonantly scattered and destroyed by Auger events (Ross, 
Fabian & Brandt 1996). Moreover, as the fraction of com- 
pletely stripped iron increases the line may be significantly 
broadened by Comptonization and so become less apparent 
in the observed spectrum (Matt, Fabian & Ross 1996). Rel- 
ativistic blurring when the reflection is from an accretion 
disc can further merge the line and edge so that neither is 
distinct in the flnal spectrum. The net result can be a weak 
broad absorption trough which starts at about 7 keV. 



Unlike the reflection spectrum from relatively neutral 
gas, which depends mainly upon the relative abundance of 
the elements and can be computed in a straightforward man- 
ner by the Monte-Carlo technique, the spectrum when the 
gas is partially ionized requires detailed numerical calcu- 
lation to obtain the temperature, ionization structure and 
resultant Comptonization (Ross 1979). We have previously 
computed examples appropriate for simple accretion discs 
around AGN and BHC, concentrating mainly on ionization 
parameters ^ = 4-KF/r^ ^ 10^ where Auger destruction can 
be important (Ross & Fabian 1993) . Since however the pre- 
cise density structure of the outer few Thomson depths of 
an accretion disc where the reflection spectrum is formed 
is unknown (the total disc thickness may be hundreds of 
Thomson depths), we consider a wider range of conditions 
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here and highlight the range at higher values of ^ ~ 10*, the 
details of which have been largely ignored in previous work. 

This range of ^ may be particularly relevant to the spec- 
tra of many BHC in the low/hard state, which show only 
a small reflection component (Ebisawa 1991; Ebisawa et al 
1994, 1996; Gierliiiski et al 1997; Zycki, Done & Smith 1998; 
Done & Zycki 1998). One popular interpretation for this is 
that the central parts of the disc out to 30 - 50 gravita- 
tional radii are missing and replaced by a hot cloud which 
gives no intrinsic reflection (Gierliriski et al 1997; Zycki et al 
1998; Poutanen 1998). What reflection is seen comes from 
the outer disc. We point out the alternative that, in a geom- 
etry where the primary X-rays are produced above a disc, 
an apparent lack of reflection could indicate high reflection. 
In principle, there is no difference in the spectrum expected 
from a single hot cloud and one with a perfect mirror bi- 
secting it. In practice, Compton reflection is not perfect, 
but departures from a uniform albedo can be small in the 
1-30 keV range where most spectral fitting takes place. We 
show here that the small departures from uniformity due to 
iron features that occur in a highly ionized disc can account 
for the small refiection signature seen, without the need to 
illuminate the outer disc. This solution may also be relevant 
to the X-ray spectra of many quasars, which also show lit- 
tle evidence for reflection yet are generally assumed to be 
powered by disc accretion (Nandra et al 1995, 1997). 

We note that values of ^ ~ 10* appear to imply an 
accretion rate close to the Eddington limit for a standard 
accretion disc (Ross & Fabian 1993 show that ^ = 7 x 10*/^, 
where / is the Eddington ratio) . However, if there is a patchy 
corona above the disc (see e.g. Stern et al 1995), rather than 
a continuous one, with a covering fraction of /a then ^ > 10* 
requires / > /a/7, which is only a per cent or so when /a ~ 
0.1. Moreover, these estimates assume a uniform density for 
the disc right to its surface. The outer few Thomson depths 
could well have a lower density and thus mean that ^ is 
underestimated. 

In this paper we present and discuss detailed, self- 
consistent computations of the temperature and ionization 
structure of slabs of gas ionized by the incident radiation and 
of the resulting reflection spectra. We assume both a con- 
stant density for the gas and a gaussian fall off. The spectra 
will be compared with X-ray observations of BHC in future 
work. 



2 COMPUTATIONS 

2.1 Method 

We consider a slab of gas whose surface is illuminated by 
radiation with a power-law spectrum of photon index F. The 
radiative transfer is calculated for the upper layer of the slab 
that is responsible for producing the reflected spectrum. In 
order to concentrate solely on the effects of the illumination, 
no radiation is taken to enter the treated layer from the 
remainder of the slab beneath it. 

The method used has been described in detail by Ross 
& Fabian (1993). The illuminating radiation is treated ana- 
lytically in a 'one-stream' approximation. The diffuse radi- 
ation that results from Compton scattering of the incident 
radiation and from emission within the gas is treated using 



the Fokker-Planck/diffusion method of Ross, Weaver & Mc- 
Cray (1978) modified for plane-parallel geometry. The local 
temperature and ionization state of the gas is found by solv- 
ing the equations for thermal and ionization equilibrium, so 
that they are consistent with the local radiation field. Hy- 
drogen and helium are assumed to be fully ionized, while 
the following ionization stages of the most abundant met- 
als are treated: C v-vil, O v-ix, Mg ix-xill, Si xi-xv, and 
Fe xvi-xxvii. 

For a given value of F, the temperature and ionization 
state of the outer portion of the slab is expected to depend 
primarily on the value of the ionization parameter, 

^-^' « 

where F is the total illuminating fiux (from 0.01 — 100 keV) 
and riH is the hydrogen number density. For uniform- density 
slabs, we vary ^ by changing the total illuminating fiux 
while keeping the hydrogen number density fixed at no = 
10^^ cm^"'. This is a typical density that might be found in 
an AGN accretion disc. At the higher density expected for a 
BHC accretion disc, the effects of three-body recombination 
and heating due to free-free absorption cause the tempera- 
ture and ionization state to depend somewhat on nn as well 
as ^ (e.g., see Ross 1979). However, this should have very 
little effect on the X-ray spectral features due to iron, which 
should remain similar to those calculated here. 



2.2 Results 

Figure 1 shows the results for a uniform slab illuminated 
by a F = 2 spectrum with an ionization parameter ^ — 10*. 
The illuminating and reflected spectra are displayed as EFe 
versus E, where E is the photon energy and Fe is the 
spectral energy flux. The slab is highly reflecting, with 
FB(out)/FE(in) > 63 per cent throughout the 2-20 keV 
spectral band. This is because the gas is highly ionized. At 
the illuminated surface, the iron is 85 per cent fully ionized, 
14per cent Fe xxvi, and 1 per cent Fe xxv. The Fe xxvi frac- 
tion peaks at 45per cent around Thomson depth tt ~ 2.5. 
Fe xxv is the dominant ion for 3.5 ^ tt ^ 7.5. Silicon and 
magnesium (not shown in Fig. 1) are fully ionized through- 
out the regions where Fe xxv-xxvil dominate. For tt ^ 10, 
iron ions have all their L-shell electrons, while magnesium 
and silicon ions have flUed K-shells. 

Despite the high ionization of the surface layers, the 
reflected spectrum shows features due to iron Ka emission 
and K-shell absorption. Most of the Ka photons emerge in a 
broad Comptonized emission feature that blends smoothly 
into the Compton-smeared absorption feature. Only a small 
fraction of the Ka photons emerge in narrow Fe xxv and 
Fe xxvi line cores at 6.7 and 7.0 keV, respectively; these 
are shown in Fig. 1 with a spectral resolution 5E / E ~ 2 per 
cent. Fe xxvi Ka photons are subject to resonance trapping. 
Many are removed from the narrow line core by an initial 
Compton scattering and then are further Comptonized as 
they diffuse outward to the surface. The Fe xxv intercom- 
bination line is not subject to resonance trapping, but most 
of these photons are produced at such great depth that they 
also suffer repeated Compton scatterings before escaping. To 
clarify the contribution of the Comptonized line. Fig. 1 also 
shows the reflected spectrum when the iron Ka emission is 
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Figure 1. Illumination of a uniform slab with F = 2 and ^ = 10*. 
Top panel: Reflected spectrum. Bottom panel; Iron ion fractions 
versus Thomson depth. 



artificially suppressed. The Ka photons form an important 
part of the reflected continuum for 4 ^ £ JS 9koV. Also 
shown in Fig. 1 is the reflected spectrum if all iron features 
are suppressed by setting the iron abundance to zero. Rela- 
tive to this smooth spectrum, the actual reflection spectrum 
is enhanced for E < 7.5 keV and depleted for E > 7.5 keV 

For soft X-rays {E < 1.5 keV in this case) not shown 
in Fig. 1, the emergent spectral flux exceeds the incident 
flux due to bremsstrahlung emission by the hot surface lay- 
ers and "inverse Compton" upscattering of even softer pho- 
tons. The total flux leaving the surface over the entire spec- 
tral range under consideration (0.01 keV < E < 100 keV) 
equals the total incident flux, as required by the condi- 
tion of thermal equilibrium. The reflected spectrum declines 
steeply above ~ 50 keV because of the sharp cutoff assumed 
at 100 keV. Extending the illuminating spectrum to higher 
energies (with an exponential cutoff, say) would raise this 
portion of the emergent spectrum via Compton downscat- 
tering of higher energy photons, but would have little effect 
on the iron ionization structure and spectral features. 

The temperature at the illuminated surface is found to 
be 1.5 X 10^ K. In the highly ionized gas there, heating due 
to Compton downscattering of hard photons dominates over 



heating due to photoabsorption. Therefore the temperature 
should not exceed the "Compton temperature," Tc, at which 
Compton heating would be balanced solely by cooling due to 
"inverse Compton" upscattering of soft photons. This con- 
dition is given by 

f--B2 



4fcTc 



ue dE ■ 



UE 



E' 



21E-' 
5moc2 



dE, 



(2) 



where ue is the spectral energy density of the radiation, 
while El and E2 are the lower and upper limits, respectively, 
of the spectrum under consideration. The right-hand side 
of this equation includes the reduction in Compton heating 
due to first-order Klein-Nishina corrections to the scattering 
cross section (see Ross 1979). Setting ue oc E^~^ for the 
illuminating radiation, this gives Tc = 1.9 x 10^ K for F = 2. 
The temperature that we find at the surface is slightly lower 
than the Compton temperature due to additional cooling by 
bremsstrahlung emission. 

This disagrees with the results of Zycki et al. (1994), 
who found a surface temperature exceeding 4 x lO'^ K for illu- 
mination with ^ = 10'' by a spectrum that was only slightly 
flatter (F = 1.9). The Monte Carlo calculation of Zycki et al. 
only treated photons in the range 0.15 keV < E < 100 keV 
and did not include the thermal emission by the gas itself. 
This leads to underestimation of the inverse Compton cool- 
ing rate and thus to the high surface temperature (and the 
steep temperature gradient) that they found. 

Figure 2 shows a series of reflection spectra under 
similar conditions with ionization parameters ranging from 
^ = 30 to ^ = 10"". The model with the highest ionization 
parameter, ^ — 10^, is an excellent reflector and exhibits 
negligibly small spectral features due to iron. This is be- 
cause the surface layer is fully ionized to great depth, with 
Fe XXVI not becoming dominant until tt ~ 8. The Compton 
reflection produces a slight steepening in the reflected spec- 
trum compared to the illumination; the reflected component 
mimics a power law with F = 2.14 in the 2-20 keV band. 
The temperature at the illuminated surface is found to be 
1.9 X lO'^ K, the full Compton temperature for a F = 2 power 
law spectrum. 

When the illuminating flux is reduced so that ^ — 
3x IC*, the Ka emission and K-shell absorption features due 
to iron begin to become apparent. Fully- ionized iron dom- 
inates for Tt ^ 5, so the emission and absorption features 
are weak and are highly broadened by Compton scattering. 
When ^ is reduced to lO**, the broad spectral features due to 
iron become more important. This model has already been 
discussed in detail. 

With ^ reduced to 3 x 10^ , less than half of the iron is 
fully ionized at the illuminated surface, and Fe xxv becomes 
dominant at tt ~ 1. Now an important narrow emission 
line due to Fe xxv can be seen in addition to a Compton- 
broadened emission feature in Fig. 2. This is primarily due 
to emission of the intercombination line following recombi- 
nation to excited states. When ^ is further reduced to lO"', 
Fe xxv dominates at the illuminated surface, and the narrow 
Fe xxv Ka line is quite strong. The tiny emission feature 
just above 8.8 keV is due to radiative recombination directly 
to the ground level of Fe xxv. Similarly, Si xiv produces 
two emission features: one at 2.01 keV due to Ka emission 
and the other just above 2.67 keV due to radiative recombi- 
nation. (See Zycki et al. 1994 for other examples of radia- 
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Figure 2. Reflection spectra for illumination of a uniform slab with F = 2 and various values of the ionization parameter. The number 
above each curve indicates the corresponding ^ value, and the dotted line to the right of that number indicates the constant value of 
EFe for that illumination. 



tive recombination emission features.) The exact strength 
of these and other low-energy spectral features could be af- 
fected by absorption by elements such as nitrogen and neon 
which we do not include in our calculations, and we do not 
calculate the features due to sulfur. 

The narrow Fe Ka emission line is suppressed in the 
models with ^ — 300 and ^ = 100. This is because ions 
in the range Fe xvil-Fe xxill dominate at the illuminated 
surface, and their Ka photons are assumed to be destroyed 
by the Auger effect during resonance trapping (see Ross & 
Fabian 1993; Zycki & Czerny 1994; Ross et al. 1996). The 
narrow line seen when ^ — 300 is due to a small amount of 
Fe XXV near the surface. On the other hand, the narrow line 
seen when f = 100 is due to Fe xvi, the least ionized species 
that we treat, which then dominates for tt ^ 0.5. Finally, 
for ^ = 30 the reflection is similar to that of a cold, neutral 
slab, and the narrow emission line at 6.4 keV is strong. 

The ionization structure in the outer layers of the illumi- 
nated slab depends on the spectral form of the illumination 
as well as on the ratio of total flux to gas density expressed 
by the parameter ^. Figure 3 shows the results with ^ — 10* 



again, but when the illuminating spectrum is a flatter power 
law with r = 1.5. Now a greater fraction of the illuminat- 
ing photons lie in the 9-20 keV range that is so important 
in producing fully photoionized iron. As a result, Fe xxvil 
dominates to greater depth (tt ^ 5.5) than when F = 2, 
and the Compton-broadened emission and absorption fea- 
tures due to iron are not as strong. 

One of the uncertainties in modelling Compton reflec- 
tion is the density structure of the illuminated gas. In the 
models presented above, the gas density has been assumed 
to be uniform with depth. This is the case, for example, in 
the standard Shakura & Sunyaev (1973) theory of accretion 
discs (also see Laor & Netzer 1989). However, this proba- 
bly is not realistic even for a bare accretion disc (e.g., see 
Shakura, Sunyaev & Zilitinkevich 1978; Shimura & Taka- 
hara 1993), and it certainly cannot be the case when the 
surface has strong external illumination. The heating of the 
outermost layers by the impinging radiation should result in 
a decrease in density there. 

In order to see the general effect of such a decrease 
in density, let us arbitrarily assume that a constant-density 
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Figure 4. Illumination with F = 2 and ^o = 10'* when there is a 
gaussian dropoff in density at the top of the slab with a Thomson 
depth TT = 5. The reflected spectrum and the incident illumina- 
tion are shown by the solid and dashed curves, respectively. Also 
shown for comparison is the reflected spectrum (dotted curve) for 
uniform density from Fig. 1 



compared to only 69 per cent fully ionized for the uniform- 
density slab shown in Fig. 1. In fact, Fe xxvil dominates all 
the way down to tt ~ 3.5. The broad, Comptonized, iron 
Ka emission and K-shell absorption features in the reflected 
spectrum are not as strong as for the uniform-density slab. 
Of particular importance is the fact that the narrow Ka line 
cores, which were already weak in the uniform-density case, 
are now almost totally suppressed. 



Figure 3. Illumination of a uniform slab with i; = 10* for a 
F = 1.5 power law. Top panel: Reflected spectrum (solid curve) 
and illumination (dashed line). Bottom panel: Iron ion fractions 
versus Thomson depth. 



slab (with nu ~ no) is topped by a "boundary layer" in 
which the density decreases with height in a gaussian man- 
ner, 



nniz) = noexp ( — — 



(3) 



Here z is the height above the base of the boundary layer, 
and h is its characteristic thickness. Since the illuminating 
radiation is expected to have an important effect down to a 
Thomson depth of a few, we let the boundary layer have a 
total Thomson depth 



Tt = 1.2noO"Tft = 5, 



(4) 



where the free electron density is assumed to be ria = 1.2nu- 
Figure 4 shows the result when such a gas is illuminated by 
a r = 2 power law with a flux that would yield an ionization 
parameter ^o = 'iivF/no — 10* for gas at the base density. 
Since the gas density is lower in the boundary layer, the 
effective ionization parameter is higher there. As a result, 
iron is 88 per cent fully ionized at Thomson depth tt = 1, 



3 DISCUSSION 

For X-ray illumination with ^ ^ IC* or higher, the features 
in the reflected spectrum due to iron Ka emission and K- 
shell absorption are weak and are smeared out by Compton 
scattering. Any narrow Fe Ka line cores are extremely weak. 
These effects are further enhanced when the illuminating 
spectrum is harder (flatter) or when there is a dropoff in 
gas density due to the heating by the external illumination. 

Such effects may come into play in the formation of X- 
ray spectra of Black Hole Candidates. In the past, Ginga 
spectra of BHCs have been interpreted as exhibiting broad 
iron absorption features ("smeared edges") with very weak 
Fe Ka lines (Ebisawa 1991; Tanaka 1991; Ebisawa et al. 
1994). This led to the suggestion that the line is suppressed 
by resonant Auger destruction in Fe xvil-xxill (Ross & 
Fabian 1993; Ueda, Ebisawa & Done 1994; Ross et al. 1996). 

However, recent spectral studies of Ginga, EXOSAT, 
and ASCA BHC in the low/hard state have found broad Fe 
Ka emission features as well as K-shell absorption features 
(Zycki, Done & Smith 1997; Zycki et al. 1998; Done & Zycki 
1998). These studies have treated the broadening of the fea- 
tures as being due to relativistic smearing by the reflect- 
ing accretion disc. The ^ values derived for the illumination 
have been very low, and the weakness of the iron features 
has been taken to imply that the disc subtends a solid an- 
gle of illuminating radiation considerably smaller than 2-k. 
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Figure 5. Illumination of a uniform slab with ^ = 10'* for a F = 
1.7 power law. Top panel: Illumination (dashed line) and reflected 
spectrum (solid curve). Bottom panel: Ratio of total spectrum 
(incident plus reflected) to best-fit power law model after taking 
into account smearing by relativistic effects for emission by an 
annulus at 7Rs around a black hole. 



(Done & Zycki conclude that the broadening of the disc Ka 
line makes it difhcult to detect with ASCA, so the narrow 
line found by Ebisawa et al. 1996 comes from the compan- 
ion star.) This has led to the conclusion that the optically 
thick accretion disc only extends inward to a few tens of 
gravitational radii. 

As an alternative, the broad and weak iron features 
could be due to illumination of the disc with ^ ^ 10'*. 
Figure 5 shows the incident and reflected spectra for illumi- 
nation with ^ — 10'' and a power law index F = 1.7. If the 
primary X-ray emission is from a corona immediately above 
the accretion disc, the reprocessor subtends a full Q — 2-k, 
and the total observed spectrum is the sum of the illumi- 
nating flux (due to the half of the radiation emitted in the 
outward direction) and the reflected flux. Fig. 5 shows the 
result after relativistic smearing if this total flux is assumed 
to originate from a narrow annulus at radius 7Rs inclined 
at 30° around a Kerr black hole. Over the 2-20 keV range, 
the best-fltting power law now has a slope of F = 1.79. 



The ratio of the observed spectrum to the best-fit power 
law model shows broad emission and absorption features 
markedly similar to those found by Done & Zycki (1998) 
for the EXOSAT spectrum of Cyg X-1. There are several 
features to note here: Comptonization of the edge means 
that it begins around 6 keV (dash-dot curve in Fig. 1) and 
now resembles more a symmetrical trough than a convential 
photelectric absorption edge; the Comptonized and smeared 
line emission fills in the lower energy part of this trough so 
that the whole feature mimics an edge at about 7 keV. 

Encouraged by the similarity bewteen our model and 
the observed spectra shown by Done & Zycki (1999), we have 
fitted a grid of our models to the brightest of the archival 
EXOSAT spectra. The models are relativistically blurred 
using the Kerr metric kernel of Laor (1990). The best fit 
over the 3-15 keV range indicates ^ — 7400 and F = 1.5 
for an iron abundance twice the Solar value (Fig. 6, upper 
and centre panels) from a disc inclined at less than 30 deg 
and extending from 1.235771 to IOOtti, with surface emissivity 
varying as (radius)"^. In the lower panel we show the ratio of 
the data to the best-fitting power law model, which strongly 
resembles the similar plot by Done & Zycki (1999) for all the 
EXOSAT d&ta.. 

It is interesting to compare the result of our calcu- 
lation with that given by the PEXRIV code (Magdziarz 
& Zdziarski 1995) in the XSPEC package for a value of 
S, = 5000, which is within the range allowed by that code 
(we correct for the different energy ranges used to define ^ 
in our codes). This is shown in Figure 7. There is clearly a 
large difference, particularly in the position of the iron edge 
(the PEXRIV code does not predict the iron line properties, 
which must be added separately). The PEXRIV code does 
not take account of the Comptonization of the features in 
the outer, most highly ionized layers. We advise against the 
use of this code for high values of ^ and for the use of our 
approach, or that of Zycki et al (1994; see also Bottcher, 
Liang & Smith 1998). 

A detailed model for an ionized accretion disc will neces- 
sarily require a range of ^ to be present, both from different 
radii, and at different distances from the ionizing source if 
the corona is patchy. The results shown in this paper will be 
useful as a guide to situations dominated by highly ionized 
matter, for example where the size of separate patches of a 
corona exceed their height above the disc. We intend to pur- 
sue a detailed comparison with observed spectra in future 
work. 
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